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<n : 

We study the deformation of a lipid membrane in response to a local pH modification. Experimentally, a basic solution is 
microinjected close to a giant unilamellar vesicle. A local deformation appears in the zone of the membrane that is closest to 
the micropipette, and relaxes when the injection is stopped. A theoretical description of this phenomenon is provided. It takes 
fully into account the spatiotemporal evolution of the concentration of hydroxide ions during and after the microinjection, as 
well as the linear dynamics of the membrane. This description applies to a local injection of any substance that reacts reversibly 
with the membrane lipids. We compare experimental data obtained in the domain of small deformations to the results of our 
I/*} 'linear description, and we obtain a good agreement between theory and experiments. In addition, we present direct experimental 
observations of the pH profile on the membrane during and after the microinjection, using pH-sensitive fluorescent lipids. 
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Introduction 

During cell life, membranes are subjected to an inhomoge- 
neous and variable environment, which can be strongly cou- 
pled to biological processes. Local pH inhomogeneities at the 
cellular scale are ubiquitous and especially important. For in- 
stance, in various types of migrating cells, an actively gener- 
ated intracellular pH gradient exists along the axis of move- 
, ment, and it appears to be essential for cell migration 1 . The 
effect of an extracellular pH gradient on cell migration and 
morphology has also been demonstrated in the case of human 
melanoma cells 2 . Local pH is also of great importance in mi- 
tochondria. The Z/+-ATP synthase enzymes that synthesize 
adenosine triphosphate (ATP), the cell's fuel, are powered by 
the local pH difference across the inner membrane of mito- 
chondria. Both the proton pumps that maintain actively this 
pH difference and the enzymes that use it are located in dy- 
namic membrane invaginations called cristae, where pH het- 
erogeneities are thought to be especially important 3 . 

In order to understand the fundamental phenomena at stake 
in the response of a biological membrane to a local pH mod- 
ification, we work on giant unilamellar vesicles (GUVs). 
These biomimetic membranes are made of amphiphilic lipid 
molecules that self-assemble into closed bilayers in water. The 
local environment of a GUV can be modified in a controlled 
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way by microinjecting a reagent close to it 4 . In particular, we 
showed previously that modifying locally the pH on a GUV 
through microinjection results into a local dynamic membrane 
deformation: the membrane shape is tightly coupled to local 
pH inhomogeneities 5 " 7 . We also developed a theoretical de- 
scription of the dynamics of a membrane subjected to a lo- 
cal modification affecting its physical properties, in the sim- 
ple case of a constant modification of the membrane involving 
only one wavelength^ 7 -. 

In this paper, we extend our theoretical description to ac- 
count fully for the time-dependent profile of the fraction of 
chemically modified lipids in the membrane. This profile re- 
sults from the diffusion of the basic solution in the water that 
surrounds the membrane. We solve analytically the corre- 
sponding diffusion problem using Green's functions, and we 
integrate numerically the linear equations describing the dy- 
namics of the membrane for the resulting time-dependent pro- 
file of the fraction of chemically modified lipids in the mem- 
brane. We compare the results of this extended theoretical 
description to experimental measurements of the height of the 
membrane deformation during and after the local injection of 
a basic solution, in the regime of small deformations. In addi- 
tion, we present a direct experimental visualization of the pH 
profile on the membrane during and after the microinjection, 
obtained using a pH-sensitive fluorescent membrane marker. 



Materials and methods 

Membrane composition and vesicles preparation 

The following lipids, from Avanti Polar Lipids, were 
used without further purification: egg yolk L-a- 
phosphatidylcholine (EYPC), brain L-a-phosphatidylserine 
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(PS), and the fluorescent lipid analog 1,2-dioleoyl-sn-glycero- 
3-phosphoethanolamine-N-(carboxyfluorescein) ammonium 
salt (DOPE-CF). All other chemicals were of highest purity 
grade: calcein and NaOH, Sigma. 

Giant unilamellar vesicles were formed by the liposome 
electroformation method 8 in a thermostated chamber. Lipo- 
some preparations for phase contrast microscopy experiments 
were made with a unique lipid mixture of EYPC and PS with 
EYPC/PS 90:10 mol/mol. For fluorescence observations, 1% 
(mol) of the fluorescent lipid analog DOPE-CF was added to 
this mixture. The particular electroformation protocol used 
in this work was the following: lipid mixture solutions were 
prepared in chloroform/diethyl ether/methanol (2:7:1) with a 
total lipid concentration of 1 mg/ml. A droplet of this lipid 
solution (1 fil) was deposited on each of the two parallel plat- 
inum wires constituting the electroformation electrodes, and 
dried under vacuum for 15 min. An AC electrical field, 10 Hz, 
0.26 Vpp, was applied to the electrodes. Water (temperature 
25°C) was added to the working chamber (avoiding agitation). 
The voltage was gradually increased (for more than two hours) 
up to 1 Vpp and maintained during 15 more minutes, before 
switching the AC field off. The GUVs were then ready for 
further use. In each preparation at least 10 GUVs of diameter 
50-80 fim were available. 

Large unilamellar vesicles (LUVs) were prepared using the 
extrusion method 9 , implemented as in Ref.— . Samples were 
prepared by dissolving and mixing the above-mentioned lipids 
in chloroform/methanol (9.4:0.6 vol/vol) to obtain the desired 
composition (EYPC/PS 90:10 mol/mol, to which 1% mol of 
DOPE-CF was then added). 

Microscopy imaging and micromanipulation 

We used a Zeiss Axiovert 200M microscope, equipped with 
a charged-coupled device camera (CoolSNAP HQ; Photo- 
metrics). The experiments were computer-controlled using 
the Metamorph software (Molecular Devices). The morpho- 
logical transformations and the dynamics of the membrane 
were followed by phase contrast and fluorescence microscopy 
(Zeiss filter set 10, Ex/Em = 475/520 nm). 

Tapered micropipettes for the local injection of NaOH were 
made from GDC-1 borosilicate capillaries (Narishige), pulled 
on a PC-10 pipette puller (Narishige). The inner diameter 
of the microcapillary used for performing the local injections 
onto a GUV was 0.3 /im. For these local injections, a mi- 
croinjection system (Eppendorf femtojet) was used. The mi- 
cropipettes were filled with basic solutions of NaOH with 
concentrations ranging from 5 to 100 mM. The injected vol- 
umes were on the order of picoliters and the injection pressure 
was 200 hPa. The positioning of the micropipettes was con- 
trolled by a high-graduation micromanipulator (MWO-202; 
Narishige). The injections were performed at different dis- 



tances from the GUV surface, taking care to avoid any contact 
with the lipid membrane. 

Steady-state fluorescence measurements 

Steady-state fluorescence measurements of LUV samples 
were carried out with a Cary Eclipse spectrofluorimeter (Var- 
ian Instruments) equipped with a thermostated cuvette holder 
(±0.1 °C). Excitation and emission slits were adjusted to 5 
nm. Fluorescence emission spectra were all recorded at 25 °C. 
All fluorescence measurements were carried out at a total lipid 
concentration of 0.2 mM. In the experiments, the pH of the 
LUV samples was gradually modified by adding aliquots of 
acid or basic solutions. The measurements were carried out 
after a few minutes of equilibration under agitation. 

Experiments 

Observation of the membrane deformation 

The chemical modification of the membrane was achieved 
by locally delivering a basic solution of NaOH close to the 
vesicle. This local increase of the pH should affect the head 
groups of the phospholipids PS and EYPC forming the exter- 
nal monolayer of the membrane—, as well as the fluorescent 
marker (when present). 

FigureQ]shows a typical microinjection experiment. We in- 
ject the basic solution during a time T = 4 s. One can see in 
Fig. Q]the vesicle before any microinjection (frame s). A 
smooth local deformation of the vesicle develops toward the 
pipette during the microinjection (first line of images). Once 
the injection is stopped, the membrane deformation relaxes 
(second line of images). This deformation is fully reversible. 
For the sake of clarity, the deformation presented in Fig. [T]is 
actually the largest this paper deals with. Indeed, we focus 
on the regime of small deformations in order to remain in the 
framework of our linear theory. In particular, it is necessary 
for our theory to be valid that the deformation height be much 
smaller than the distance between the membrane and the mi- 
cropipette. 

Observation of the pH profile on the membrane 

We present here a direct experimental visualization of the pH 
profile on the membrane during and after the local microin- 
jection of NaOH. For this, we use a pH-sensitive fluorescent 
membrane marker (DOPE-CF). The fluorescein group is a 
weak acid whose conjugate base has a strong fluorescence, 
and which is attached to phosphatidylethanolamine lipids. 

First of all, the dependence of the fluorescence on the pH 
was verified on a solution of LUVs (see Fig. |2j. The depen- 
dence of the mean intensity of the LUV solution on the pH 
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Fig. 1 Typical microinjection experiment, lasting T = 4 s, observed using phase contrast microscopy. A local modulation of the pH on the 
vesicle membrane induces a smooth deformation of the vesicle (frames 0.33 to 4.08 s). The deformation is completely reversible when the 
NaOH delivery is stopped (frames 4.21 s to the end). 



is well-described by the sigmoidal shape characteristic of an 
acid-base titration (see Fig. The steep intensity rise is ob- 
served around pH 7.5, which makes this marker adequate to 
investigate pH increases starting from a neutral pH. 




Fig. 2 Mean (fluorescence) intensity as a function of the pH in a 
solution of LUVs whose membranes contain a pH-sensitive 
fluorescent marker (DOPE-CF). Dots: experimental data. Plain line: 
sigmoidal fit. 

Using this fluorescent marker within a preparation of 
GUVs, it is possible to observe the pH profile on the mem- 
brane during and after the local microinjection of NaOH. Typ- 



ical results are presented in Fig. [3] The pH profile on the 
membrane is visualized directly together with the deforma- 
tion of the membrane. The vesicle deforms progressively in 
response to the local pH increase (frames from to 3.6 s). 
The increase of the intensity in front of the micropipette and 
the lateral spreading of the bright zone during the injection 
are visible in the three first images. This illustrates the lo- 
cal pH increase on the membrane in front of the micropipette, 
and then on the sides, as the HO ions diffuse from the mi- 
cropipette tip towards the membrane. Later on, photobleach- 
ing occurs, and its effect is visible faster in the zone where the 
onset of fluorescence occurred sooner, i.e., just in front of the 
micropipette. The deformation relaxes fully when the NaOH 
injection is stopped (frames from 4.4 s to the end), and the 
fluorescence decreases at the same time, both due to the diffu- 
sion of the injected basic solution in the water surrounding the 
membrane after the end of the injection, and to photobleach- 
ing. 

The membrane shape profiles and the membrane intensity 
profiles have been extracted from the images in Fig. [3] using 
the software Image J supplemented by our own plugins. The 
corresponding results are shown in Fig. |4]for the frames cor- 
responding to times 0, 0.8, 1.2 and 2.8 s of Fig. [3] One can see 
first the local increase of intensity, and thus of pH, in front of 
the micropipette, as the membrane deforms (red curves), and 
the subsequent spreading of the high-pH zone while the pH 
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Fig. 3 Typical microinjection experiment, lasting T = 4 s, observed using fluorescence microscopy. The white arrow represents the 
micropipette tip. The pH profile on the membrane is visualized directly together with the deformation of the membrane. The vesicle deforms 
progressively in response to the local pH increase (frames from to 3.6 s). The increase of the intensity illustrates the pH increase. The 
deformation relaxes fully when the NaOH injection is stopped (frames from 4.4 s to the end). The initial fluorescence of the membrane at 
t = 0, i.e., at pH 7, is weak but visible (see Fig.(2]l. The global darkening of the images with time is due to photobleaching. 



increase continues (blue curves). In the last (green) curves, 
the effect of photobleaching is visible, as the central intensity 
decreases while the injection still continues and the membrane 
keeps deforming. 




x((jm) 



Fig. 4 (a) Membrane shape profiles during the experiment presented 
in Fig.|3](black: t=0 s, red: t=0.8 s, blue: t=1.2 s, green: t=2.8 s). 
Insert: zoom of the central zone, (b) Intensity profiles on the 
membrane at the times 0.8, 1.2 and 2.8 s (same colors as in (a)). 



Theoretical description 

Membrane free energy and force density 

In our experiments, only the external monolayer is directly 
affected by the chemical modification. Indeed, the perme- 
ation coefficient of hydroxide ions through a lipid membrane 
is Pq H - ~ 10~ 5 cm/s, which yields a negligible pH increase 
inside the vesicle on the timescale of our experiment s U ' 12 . 
Thus, to study the membrane deformation, it is necessary to 
take fully into account the bilayer structure of the membrane. 
The free energy per unit area of a bilayer can be written as 
the sum of the free-energy densities of the two monolayers, 
which will be noted /+ and f~ . We shall use a description of 
the membrane based on a local version of the area-difference 
elasticity model 1314 . The local state of each monolayer is de- 
scribed by two variables: the total curvature c defined on the 
membrane midlayer (c is thus common to both monolayers), 
and the scaled density = (p ± — po) / po defined on the mid- 
layer, where po is a reference density. The free energy per unit 
area of a monolayer reads 

= ? + ^±^c + k -^±ec)\ (1) 

where Ob 1S the membrane tension and k is its bending rigidity, 
while k denotes the stretching modulus of a monolayer, and e 
stands for the distance between the neutral surface 15 of each 



monolayer and the midsurface of the bilayer (see Fig. [5j. All 
these constants are identical for the two monolayers because 
they are assumed to be identical before the injection. The two 
monolayers also have opposite spontaneous curvatures, noted 
=Fco, since their lipids are oriented in opposite directions. We 
choose the sign convention for the curvature in such a way 
that a spherical vesicle has c < 0. With this convention, the 
monolayer denoted "+" in Eq. (fl~|i is the external monolayer. 

The expression (Q]i of the free-energy density of a mono- 
layer is a general second-order expansion around a reference 
state characterized by a fiat shape (c = 0) and a uniform den- 
sity p ± — pcr^- It is valid for small deformations around this 
reference state: r — &{e) and ec = &{e), where £ is a small 
nondimensional parameter characterizing the deformation. 

Let us call the mass fraction of the lipids of the exter- 
nal monolayer that are chemically modified by the local pH 
increase in the experiment (see Fig. [5}- Given the small con- 
centrations used and the micropipette-membrane distance, the 
pH on the membrane remains well below the effective pK a of 
the lipid head groups (which are of order 10 or higher 16 ' 17 ). 
Hence, we assume that <j> = &(e). In addition, as justified 
above, we consider that the inner monolayer is not affected by 
the pH increase outside the vesicle. Including the small vari- 
able in addition to r + and ec, we can write to second order 
in e: U& 

/+ = ^ + ( j 1 (/) + y0 2 + d(l+r + )(/)ln(/) + | C 2 

+ |(c + c </>)c + ^(r + +ec) 2 . (2) 

Note that the term in <p In in Eq. (ffjl corresponds to mixing 
entropy 14 . The constants <7i, Oi, & an d cq describe the re- 
sponse of the monolayer to (j). 

To first order, the chemical modification of the external 
monolayer (i.e., the presence of a nonzero 0) causes both a 
change of the equilibrium density and a change of the spon- 
taneous curvature of this monolayer 7 . On the neutral surface, 
these two effects are decoupled 15 : let us thus discuss them 
on this surface. For a homogeneous monolayer with constant 
mass, minimizing the free energy per unit mass f + /p + with 
respect to r + and c yields a spontaneous curvature changed by 
the amount 8cq = co<p to first order, with cq = co + 2a\e/ K. 
This minimization also gives a scaled equilibrium density on 
the neutral surface changed by the amount 5r+ = G\<p/k to 
first order. 

The force density in each monolayer of a bilayer with lipid 
density and composition inhomogeneities described by the 
free-energy densities ([TJ-© has been derived in Ref.— . As 
we shall focus on small deformations of the membrane with 
respect to the plane shape, it is convenient to describe the 
membrane in the Monge gauge, i.e., by its height z = h(x,y) 
with respect to a reference plane, x and y being Cartesian co- 
ordinates in the reference plane. Then, c = V 2 /z + ^(e 3 ) for 




Fig. 5 Sketch of the situation described (not to scale). Due to the 
injection of a basic solution from the micropipette (red) standing at 
zo above the membrane, some lipids are chemically modified (dark 
blue) in the external monolayer. The mass fraction of these modified 
lipids is denoted by (j). The membrane deforms because of this local 
chemical modification: the shape of its midlayer S" is described by 
h{r). The variables c = V 2 /j and used in our theory are defined 
on y, which is at a distance e from the neutral surface Jf of 
monolayer ±. 

small deformations such that d[h = and didjh — 
where i,j G The force density in the membrane then 

reads to first order in e 

pf(x,y) = -kd i {r+ + eV 2 h-^ij , (3) 

py(x,y) = -kd i (r--eV 2 h) , (4) 

p z (x,y) = (JoV 2 /i - kV 4 h - keV 2 r a - - a^ej V 2 , 

(5) 

where pf is the force density in monolayer "±" acting in a 
direction i tangential to the membrane, while p z = pf + p7 is 
the total normal force density in the membrane. In these for- 
mulas, we have introduced the antisymmetric scaled density 
r a = r + — r~ , and the constant it = K + 2ke 2 . 

Eq. (O shows that the equilibrium density change 5r+ = 
<5\§jk yields a tangential force density proportional to V0. 
Besides, Eq. (01 shows that both the spontaneous curvature 
change 5cq = c~o0 and the equilibrium density change 5r+ 
yield a normal force density proportional to V 2 (see also 
Refs .-L!i). Hence, both effects can drive a deformation of the 
membrane, but only the equilibrium density change can in- 
duce a significant tangential lipid flow. 

Membrane dynamics 

The linear dynamical equations describing the joint evolution 
of the membrane shape h and of the lipid densities were 
first developed in Refs. 13 > 18 . We have adapted them to our 
study of the membrane deformation caused by a local chem- 
ical modification, first in the case where co = 6 , and then in 
the general case where both the equilibrium density and the 
spontaneous curvature of the external monolayer are changed 
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due to the chemical modification 7 . In the present work, we 
use the same dynamical equations as in Ref. 7 , but we now take 
fully into account the spatiotemporal profile of during and 
after the microinjection, while Refs. 6,7 focused on the simple 
case of a constant, one-wavelength, 0. 

Let us recall briefly the derivation of the linear dynamical 
equations describing the membrane dynamics 7 - 13 . One uses a 
normal force balance for the bilayer, which involves the nor- 
mal elastic force density p z in Eq. (0 and the normal viscous 
stress exerted on the membrane by the water that surrounds 
it. Besides, as each monolayer is a two-dimensional fluid, one 
writes down generalized Stokes equations involving the mono- 
layer tangential elastic force densities in Eqs. (0 or (@), the 
two-dimensional viscous stress associated with the lipid flow, 
the tangential viscous stress exerted by the surrounding wa- 
ter, and the force density associated with the intermonolayer 
friction. The latter reads ^fb(v + — v~), where b is the inter- 
monolayer friction coefficient 18 , while v ± (x,y,t) is the two- 
dimensional velocity of the lipids within monolayer ±. Since 
= 0(e), and since the lipid flow is induced by the chemi- 
cal modification, one has |v | = ff{e). Finally, one uses mass 
conservation in each monolayer. Considering that each mono- 
layer has a fixed total mass is justified here since the timescales 
of our experiments (about 10 seconds) are much shorter than 
the flip-flop characteristic time, which is assumed not to be 
significantly modified by the local chemical modification we 
study. 

The resulting equations can be expressed conveniently in 
terms of the in-plane Fourier transforms of the fields h, r and 
0, denoted with a hat, e.g., 



describes the dynamical response of the membrane 6 i 7 - 13 ; 



h(q,t) 



drh(r,t)e 



-,qr 



(6) 



where r = (x,y) and q — (q x ,q y ) are two-dimensional vectors. 
The geometry of the problem we are studying features a cylin- 
drical symmetry, the axis of revolution being that of the mi- 
cropipette. Let us take the origin of the in-plane coordinates 
(i.e., x and y) on the axis of the micropipette. Then, the fields 



h, and only depend on r 
transforms only depend on q = 



y 'x 2 + y 2 , and their Fourier 



Combining the above-mentioned dynamical equations 
yields the following first-order linear differential equation on 
the two-dimensional variable X(q,t) = (qh(q,t), f a (q,t))'-— 



-£(q,t)+M(q)X(q,t)=Y(q,t), 



(7) 



where t is time and q is the norm of the two-dimensional 
wavevector. In Eq. 0, we have introduced the matrix which 



M(q) 



kq 3 



V 



4tj 

keq 3 



keq 2 \ 
4tj 



kq 2 
~2b 



(8) 



/ 



where rj is the viscosity of water. We have assumed that 
Tl2q 2 <C b and r\q<^b, where r\2 is the two-dimensional vis- 
cosity of the lipids, which is very well verified for all the 
wavevectors with significant weight in 0. Indeed, in the exper- 
imental situation we wish to describe, the modified lipid mass 
fraction has a smooth profile with a characteristic width 
larger than 1 /xm, and typical values of the dynamical param- 
eters, used throughout, are 772 = 10~ 9 Js/m 2 , b = 10 9 Js/m 4 
(see Refs. 19 i 20 ), and r\ = 10~ 3 Js/m 3 for water. Besides, the 
forcing term in Eq. (0 reads 7 : 



Y(q,t) = 



( Kcpq 2 
877 



(9) 



Eqs. (00 show that the antisymmetric density is coupled 
to the height of the membrane deformation. Physically, it is 
the symmetry breaking between the monolayers that causes 
the deformation of the membrane. In the present case, this 
symmetry breaking is caused by the chemical modification of 
certain membrane lipids in the external monolayer, i.e., to the 
presence of 0. And indeed, Eq. (0 shows that the forcing term 
in Eq. is proportional to (j)(q,t). 

The present theoretical description is general and applies to 
any local chemical modification of the membrane. Indeed, we 
have not made any specific assumption on (j), apart that it is 
small. What changes with the nature of the reagent injected is 
the value of the constants (7i and co, which describe the linear 
response of the membrane to the chemical modification. 

Calculation of (j>(q, t) 

The time evolution of the membrane deformation during and 
after the microinjection can be determined by solving Eq. (0. 
To this end, we first need to determine (j)(q,t), which is in- 
volved in the forcing term Y(q,t) (see Eqs. (0-0). The 
acid-base (and complexation) reactions that occur between 
the lipids and the injected hydroxide ions are reversible and 
diffusion-controlled—. Hence, the local mass fraction 0(r,/) 
of chemically modified lipids is determined by an instanta- 
neous equilibrium with the hydroxide ions that are above the 
membrane. Let us denote by c (r, z, t) the local concentration in 
hydroxide ions in the water that surrounds the vesicle, and let 
us assume that the unperturbed flat membrane is in the plane 



z = 0, and that the micropipette stands at zo > above it (see 
Fig. [5j. In the experimental conditions, the pH on the mem- 
brane remains well below the effective pKa of the lipids. Thus, 
the mass fraction of chemically modified lipids after the re- 
action is proportional to the concentration of HO ions just 
above the membrane: <f>(r,t) °c c(r,z = 0,f). In this linear 
regime, calculating <j>(r,t) amounts to evaluating c(r,z = 0,t). 

In our experiments, the NaOH solution is injected from a 
micropipette of inner diameter d = 0.3 jJ.m and length L ~ 2 
cm, with an injection pressure AP = 200 hPa. Hence, we can 
estimate the average velocity vo of the NaOH solution when it 
just gets out the pipette, treating the flow in the micropipette 
as a Poiseuille flow: vo = APd 2 / \2>2r]L) ~ 2/im.s -1 . Thus, 
given the small lengthscales and velocities involved in the mi- 
croinjection, the Peclet number is very small: Pe = vqzo /D s» 
10~ 2 , where the order of magnitude of zo is 10/im, while D is 
the diffusion coefficient of sodium hydroxide in water, which 
is given by D = 2/(l/D OH - + l/£> Na +) = 2125jim 2 /s at in- 
finite dilution and at 20°C 22 . Hence, once the NaOH solution 
is out of the pipette, its dynamics is dominated by diffusion. 
This fact can also be verified experimentally by injecting a 
fluorescent substance in the conditions under which the basic 
solution is injected in our experiments. The observed quasi- 
spherical fluorescent "cloud", presented in Fig. |6] illustrates 
well the dominant effect of diffusion. 

Thus, we can consider that the NaOH solution simply dif- 
fuses from the micropipette tip (r = 0,zo). Besides, since the 
membrane is a surface and <C 1, it is possible to neglect the 
hydroxide ions that react with the membrane when calculating 
c. Given these simplifications, c can be obtained by solving 
the (three-dimensional) diffusion equation 



d,c-DV 2 c = S, 



with the source term 



S(r, z, t ) = So 8 (r) <5 (z - zo) 1 [o,r] (?) , 



(10) 



(11) 



where 1 p j-j is the indicator function of the interval [0,T]. This 
source term corresponds to a constant injection flow from the 
point (r = 0,zo) for < t < T, where we have chosen the 
injection start as our time origin, and we have called T the time 
when the injection stops. In addition, the membrane imposes 
a Neumann boundary condition to c: d z c (r,z = h(r,t),t) = 0. 
Since we are working at first order in e, we simply use 



d z c(r,z = 0,t)=0. 



(12) 



The solution to this diffusion equation with the boundary 
condition Eq. ( TTZb can be written as: 



c(r,z,t) = fdt' [ dr' f + dz' S fr' ,z' ,?') G (\r—r > \,z,£,t — t f ) 
Jo Jr 2 Jo 

(13) 



/■min(/, r) 

S J dt'G(r,z,z Q ,t-t') , 



where G is the causal Green function of the diffusion equation 
which verifies the Neumann boundary condition Eq. (fT2l . It 
can be expressed simply using the method of images 23 : 



G(r,z,z',t)^G„ (r,z-z',t)+G„ (r,z + z',t) 



(14) 



where we have introduced the infinite-volume causal Green 
function of the diffusion equation 



G^{r,z,t) 



e(t) 



exp 



r 2 + z 2 
ADt 



(15) 



{AnDtf 2 

where 9 denotes Heaviside's function. Using Eqs. (fT~3l IT4l 
[TBI , it is straightforward to obtain an analytical expression for 
c(r,z,t), and for its in-plane Fourier transform c(q,z,t). Since 
<p(q,t) is proportional to c (q,z = 0,f), we also obtain its ana- 
lytical expression, which reads 



(q,t) = f (q,t) - 0{t-T) (q,t- T) , (16) 
cosh(gzo) sinh(^zo) 



with 



2VDt 



qzo 



qzo 



(17) 

where erf denotes the error function. The value of the pro- 
portionality constant in this expression is not crucial for our 
study since all our dynamical equations are linear: it only af- 
fects the deformation and the antisymmetric scaled density by 
a multiplicative constant. 

Resolution of the dynamical equations 

Now that we have determined <j)(q,t), which is involved in the 
forcing term Y(q,t) of the dynamical equations (see Eqs. ([7]- 
[8)), the time evolution of the membrane deformation during 
and after the microinjection can be determined by solving the 
differential equation Eq. ©. This can be done thanks to the 
method of variation of parameters. The square matrix M de- 
fined in Eq. (0 has two real positive and distinct eigenvalues 
for all q > 0. Let us call these eigenvalues ji an d 72, and 
let us introduce the associated eigenvectors V\ = (vi,wi) and 
V2 = (v2,W2). For the initial condition X(q,t = 0) = (0,0), 
corresponding to a non-perturbed membrane (i.e., flat and with 
identical density in the two monolayers), we can write: 

qh(q,t) = J ds [vj e^'AfV) + v 2 e~ r2 'fi(s)] , (18) 

where A(t) and B(t) are the solutions of the linear system 

ViA(t)e- r>t +V 2 B(t)e-K' =Y{q,t). (19) 

Then, in order to obtain the membrane deformation profile 
at time t, we perform an inverse Fourier transform: 



h(r,t) 



1 

2% Jo 



dq [qh(q,t)]j {qr), 



(20) 
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Fig. 6 Injection of a water solution containing a fluorescent marker (calcein 96 mM) from a micropipette in the conditions under which the 
basic solution is injected in our experiments: injection pressure AP = 200 hPa, microcapillary inner diameter 0.3 /xm, injection duration 
r = 4s. 



where we have used the cylindrical symmetry of the problem, 
introducing the Bessel function of the first kind and of zero 
order Jo- Thus, using Eq. ( fTSI l. we finally obtain 

*(r,0 = T- l + dqMqr) f ds [vi e^'A(s) + V2<T*' B(s)] . 

27T JO JO 

(21) 

Hence, we can obtain the spatiotemporal evolution of the 
membrane deformation during and after the microinjection by 
carrying out the integrals in Eq. (|2H numerically. 

Comparison between theory and experiments 

In the previous section, we presented a theoretical description 
of the membrane deformation in response to the microinjec- 
tion of a basic solution close to a membrane. In order to com- 
pare the predictions of this description to experimental results, 
we measured the height H(t) = h(r — 0,f) of the membrane 
deformation in front of the micropipette during the microin- 
jection experiments described in the experimental section (see, 
e.g., Fig. [l}. We carried out several microinjection experi- 
ments on different GUVs, all with an injection lasting T = 4 s. 
For each GUV, several such experiments were conducted, with 
various distances zo between the micropipette and the mem- 
brane. 
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Given that our theory is linear, we expect it to be valid in 
the regime where (j) <C 1 and H <C Zo- Besides, as the unper- 
turbed membrane is considered flat in our theory, its domain 
of validity is restricted to zo <C R, where R is the radius of 
the GUV. Hence, we strived to remain in these conditions, 
and we shall present here only the experimental results that 
match these conditions best. In practice, the radii of our largest 
GUVs were of order 60 to 80 jim. We thus focused on values 
of zo in the range 10 to 30 jim. Besides, we adjusted the HO 
concentration for the various zo in order to obtain small but 
observable deformations, of order 1 to 5 |im. As long as we 
remain in the linear regime, the absolute value of the concen- 
tration of the injected solution only affects the deformation by 
a global proportionality constant. This was checked experi- 
mentally, albeit in a rough fashion given the uncertainty on 
the radii of the different micropipettes used. 

In order to compare the experimental data on H{t) during 
and after the microinjection to the solutions of our theoretical 
equations, we normalize our experimental data on //(f) by the 
value of H(t = 4s), corresponding to the end of the injection, 
for each experiment. This eliminates the effect of the unknown 
proportionality constant in our theoretical expression of (see 
Eqs. (fToT - fTTI i), as well as the experimental effect of the differ- 
ent concentrations and of the slightly different micropipette 
diameters. 



The membrane deformation H(t) = h(r = 0,f) predicted 
theoretically using Eq. ( f2TT > together with Eq. ([19) and 
Eqs. (fToT - TTTI i was calculated numerically. This was done in 
the case of four-second injections, with the values of zo corre- 
sponding to the different experiments, and taking typical val- 
ues of the membrane constitutive constants: k = 10 19 J, k = 
0.1 N/m, b = 10 9 J.s.m~ 4 ande = 1 nm, and r\ = 10~ 3 J.s.rrT 3 
for the viscosity of water. 

In order to solve the dynamical equations, it was also neces- 
sary to assign a value to the parameter a = —kecol<5\, which 
quantifies the importance of the change of the spontaneous 
curvature relative to the change of the equilibrium density of 
the external monolayer as a result of the chemical modifica- 
tion. This parameter cannot be determined from an analysis of 
static and global modifications of the environment of the vesi- 
cle. Indeed, the area-difference elasticity model predicts that 
the equilibrium shape of a vesicle is determined by a combined 
quantity which involves both the equilibrium density and the 
spontaneous curvatur e 17 ' 24 . A rough microscopic lipid model 
based on geometry yields a ss 1—: the effect of the change of 
the spontaneous curvature and that of the equilibrium density 
should have the same order of magnitude. In the absence of 
any experimental measurement of a, we took a = 1 in our cal- 
culations. We also checked that the agreement between theory 
and experiment was not as good for a =0.1 and a = 10 as it 
is for a = 1 . 

Figs.|7]and[8]show the results obtained for several values of 
Zq, on two different vesicles. As these vesicles have the same 
lipid composition, their membranes share the same constitu- 
tive constants, as well as the same constants (7i and cq that 
describe their response to (p. On the contrary, the membrane 
tension Ob is highly variable among vesicles. We did not mea- 
sure this tension during the experiments, but since the vesicles 
are flaccid, ay should be in the range 10~ 8 - 10~ 6 N/m. It 
was visible that the vesicle corresponding to Fig. |7]was more 
flaccid than that corresponding to Fig. [8] The parameter Go 
was adjusted in these two cases, with all the other parameters 
kept constant at the above-mentioned values. More precisely, 
we integrated numerically our dynamical equations assuming 
various values of the tension Gq. For each of the two vesicles, 
we calculated the total chi-square between these numerical re- 
sults and the three experimental data sets that respected best 
the hypotheses of our theory, i.e., small deformation and sta- 
ble injection pressure (see Figs. |7ja) and|HJa)). For the vesi- 
cle corresponding to Fig. [7] the best match between theory 
and experiments, i.e., the lowest chi-square, was obtained for 
<To = (1.5 ±0.5) x 10~ 8 N/m, and for the vesicle correspond- 
ing to Fig. [8] it was obtained for Ob = (5 ± 1) x 10~ 7 N/m. 
These values are in the expected range, and the data corre- 
sponding to the most flaccid vesicle is best matched by the 
lowest tension, which is satisfactory. 

The experimental data presented in Figs.|7]and|8]is slightly 
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Fig. 7 Normalized height H(t) /H(t = 4 s) of the membrane 
deformation, during and after microinjections lasting T = 4 s. The 
microinjections were carried out on the same GUV at different 
distances zo, corresponding to the different colors, (a) Comparison 
between experimental data and theoretical calculations for three 
experiments. Dots: raw experimental data (one data point was taken 
every 57 ms). Plain lines: numerical integration of our dynamical 
equations with Gq = 1.5 X 1CP 8 N/m. (b) Full set of experimental 
data; a moving average over 4 successive points was performed to 
reduce the noise, (c) Full set of numerical data for values of zo 
corresponding to those of the experiments. 



noisy. This is due to the fact that we have focused on small 
deformations, of 1 to 5 fim, in order to remain in the domain 
of application of our linear theoretical description. For such 
small deformations, all sources of noise (e.g., vibrations) be- 
come important, and the pixel size also becomes limiting. Our 
image treatment plugins were equipped with subpixel resolu- 
tion in order to improve this point (see, e.g., Ref. 25 ). Besides, 
it can be noted that the injection phase is more noisy than the 
relaxation phase, especially on Fig. [8] This is due to fluctua- 
tions of the injection pressure, which seem to occur mostly at 
the beginning of the injection phase. In particular, the exces- 
sive overshoots observed at the beginning of the injection on 
Fig. [8] are very likely due to this artifact. 

Figs. |7] and [8] show a good agreement between our experi- 
mental data and the results of our theoretical description. In 
particular, our theory predicts the right timescales of defor- 
mation and of relaxation, and also the right variation of these 
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Fig. 8 Similar data as on Fig. [7] but for another, less flaccid, vesicle. 
Here, the numerical integration of our dynamical equations was 
carried out for (To = 5 x 1CT 7 N/m. 



timescales with the distance zo between the membrane and the 
micropipette. The increase of the timescales with zo comes 
from two different factors. First, the diffusion of the HO 
ions takes longer if zo is larger. Second, when zo increases, 
the width of the modified membrane zone that deforms in- 
creases, so that smaller wavevectors have a higher importance 
in h(q,t). As the relaxation timescales of the membrane, 
which correspond to the inverse of the eigenvalues of M, all 
increase when q decreases (see Eq. (0), this yields longer 
timescales. 

Besides, the timescales involved are shorter if the vesicle 
is more tense, for instance they are shorter in Fig. |8]than in 
Fig. [7] This can be understood as follows: one of the relax- 
ation timescales of the membrane, which corresponds to the 
inverse of one of the eigenvalues of M (see Eq. ([8])), can be 
approximated by 4 77 / Ooq for the wavevectors q with largest 
weight in h(q). This timescale decreases when Co increases. 
More qualitatively, a tense membrane will tend to relax faster 
once it has been deformed. 

While we have focused our discussion on the height H(t ) = 
h(r = 0,f) of the membrane deformation in front of the mi- 
cropipette, the full deformation profiles h(r,t) are available 
both experimentally and theoretically in our work. We ob- 
serve qualitative agreement between theory and experiment 



regarding the spatial profile of h (and 0), but the effect of vesi- 
cle curvature makes it difficult to push further the quantitative 
analysis away from the micropipette axis. 

Conclusion 

We have studied experimentally and theoretically the defor- 
mation of a biomimetic lipid membrane in response to a lo- 
cal pH increase obtained by microinjecting a basic solution 
close to the membrane. Experimentally, we have measured 
the deformation height during and after the injection, and we 
have observed directly the pH profile on the membrane at the 
same time, using a pH-sensitive fluorescent membrane marker. 
Theoretically, our description of the phenomenon takes into 
account the linear dynamics of the membrane, and it also 
fully accounts for the time-dependent profile of the fraction of 
chemically modified lipids in the membrane. This profile re- 
sults from the diffusion of the basic solution in the water that 
surrounds the membrane during and after the microinjection. 
We have compared experimental data regarding the height of 
the deformation to the results of our theoretical description, 
in the regime of small deformations, and we have obtained a 
good agreement between theory and experiments. 

Experimentally, it would be interesting to measure the vesi- 
cle tension through a micropipette 26 at the same time as the 
microinjection is performed. It would perhaps become pos- 
sible to adjust the intermonolayer friction coefficient b (and 
hence to measure it) if the tension was known precisely. How- 
ever, since high tensions yield shorter timescales and smaller 
deformations, it would be necessary to control precisely very 
small tensions through the micropipette. 

From a theoretical point of view, our description is general 
and applies to a local injection of any substance that reacts 
reversibly with the membrane lipids. It could be improved 
by taking into account the curvature of the vesicle instead of 
taking a flat membrane as a reference state. Besides, it would 
be interesting to include nonlinear effects in order to describe 
larger, more dramatic, deformations, such as tubulation 6 . 

The study of the response of a membrane to a local modifi- 
cation of its environment is a promising field. From the point 
of view of membrane physics, studying the spatiotemporal re- 
sponse of a membrane to a local modification can give access 
to membrane properties that were not accessible before. For 
instance, the ratio of the spontaneous curvature change to the 
equilibrium density change caused by a chemical modification 
cannot be determined from an analysis of static and global 
modifications^. In the theoretical work Ref. 27 , we present a 
way of determining this ratio from the study of the dynami- 
cal response of a membrane to a continuous local injection. 
Hence, an interesting perspective would be to study contin- 
uous injections experimentally, which could yield an experi- 
mental measurement of this ratio. However, the injection pres- 
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sure would have to be very stable for this to be done satisfac- 
torily. More generally, we hope that studying the response of 
a biomimetic membrane to a local modification of its environ- 
ment will help to shed light onto the relation between cellular 
phenomena and small-scale environment changes. 
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